.

library(tidyverse)
library(dplyr)
library(forcats)
library(gmodels)
library(haven)

###Read data

#Put the complete name of the countries. 
Disc <- Disc |> mutate(country = case_when(
  isocntry== "BE" ~ "Belgium",
  isocntry== "DK" ~ "Denmark",
  isocntry== "GR" ~ "Greece",
  isocntry== "ES" ~ "Spain",
  isocntry== "FI" ~"Finland",
  isocntry== "FR" ~ "France",
  isocntry== "IE" ~ "Ireland",
  isocntry== "IT" ~ "Italy",
  isocntry== "LU" ~ "Luxembourg",
  isocntry== "NL" ~ "Netherlands",
  isocntry== "AT" ~ "Austria",
  isocntry== "PT" ~ "Portugal",
  isocntry== "SE" ~ "Sweden",
  isocntry== "DE-W" ~ "Germany (specifically the state of North Rhine-Westphalia)",
  isocntry== "DE-E" ~ "Germany (specifically the state of Berlin)",
  isocntry== "GB" ~ "United Kingdom",
  isocntry== "BG" ~ "Bulgaria",
  isocntry== "CY" ~ "Cyprus",
  isocntry== "CZ" ~ "Czech Republic",
  isocntry== "EE" ~ "Estonia",
  isocntry== "HU" ~ "Hungary",
  isocntry== "LV" ~ "Latvia",
  isocntry== "LT" ~ "Lithuania",
  isocntry== "MT" ~ "Malta",
  isocntry== "PL" ~ "Poland",
  isocntry== "RO" ~ "Romania",
  isocntry== "SK" ~ "Slovakia",
  isocntry== "SI" ~ "Slovenia",
  isocntry== "HR" ~ "Croatia",
  TRUE ~ NA_character_))

view(Disc)

Convert numeric answers from qc19 to character answers

Disc <- Disc %>% mutate(change_docs = case_when(
  qc19 == 1 ~ "Yes",
  qc19 == 2 ~ "No", 
  qc19 == 3 ~ "DK", 
  TRUE ~ NA_character_)) 

Count the qc19 answers by country

library(dplyr)

#sum counts 
opinions_by_country <- Disc %>%
  group_by(country, change_docs) %>%
  summarise(count = n()) %>%
  pivot_wider(names_from = change_docs, values_from = count, values_fill = 0)

opinions_by_country

#percentage counts 
Disc %>%
  group_by(country, change_docs) %>%
  summarise(count = n()) %>%
  ungroup() %>%
  group_by(country) %>%
  mutate(total_count = sum(count)) %>%
  mutate(percentage = (count / total_count) * 100) %>%
  select(country, change_docs, percentage) %>%
  pivot_wider(names_from = change_docs, values_from = percentage, values_fill = 0)
#Top 5 of countries that think that transgender persons should be able to change their civil documents to match their inner gender identity
country_with_most_yes <- opinions_by_country %>%
  filter(Yes > 0) %>%
  arrange(desc(Yes)) %>%
  head(5)
country_with_most_yes
#Top 5 of countries that do not think that transgender persons should be able to change their civil documents to match their inner gender identity
country_with_most_no <- opinions_by_country %>%
  filter(No > 0) %>%
  arrange(desc(No)) %>%
  head(5)
country_with_most_no
#Top 5 of countries that do not think that transgender persons should be able to change their civil documents to match their inner gender identity
country_with_most_DK <- opinions_by_country %>%
  filter(DK > 0) %>%
  arrange(desc(DK)) %>%
  head(5)
country_with_most_DK
LS0tCnRpdGxlOiAiRmluYWwgQXNzaWdubWVudCIKYXV0aG9yOiAiVmFsZXJpYSBDb250cmVyYXMiCmRhdGU6ICJgciBTeXMuRGF0ZSgpYCIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICB0b2M6IHllcwogICAgZGZfcHJpbnQ6IHBhZ2VkCiAgaHRtbF9ub3RlYm9vazoKICAgIHRvYzogZmFsc2UKICAgIGRmX3ByaW50OiBwYWdlZAotLS0KCi4KCmBgYHs9aHRtbH0KPHN0eWxlPgpib2R5IHsKdGV4dC1hbGlnbjoganVzdGlmeX0KPC9zdHlsZT4KYGBgCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobz1ULCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFLCBrbml0ci5wdXJsLmlubGluZSA9IFRSVUUgKQpgYGAKCmBgYHtyfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShkcGx5cikKbGlicmFyeShmb3JjYXRzKQpsaWJyYXJ5KGdtb2RlbHMpCmxpYnJhcnkoaGF2ZW4pCmBgYAoKIyMjUmVhZCBkYXRhIApgYGB7cn0KRGlzYyA8LSByZWFkX2R0YSgiRGF0YS9aQTc1NzUuZHRhIikKaGVhZChEaXNjKQoKRGlzYyB8PiBzZWxlY3QoaXNvY250cnksY291bnRyeSxxYzE5KQp1bmlxdWUoRGlzYyRpc29jbnRyeSkKCmBgYAoKCmBgYHtyfQojUHV0IHRoZSBjb21wbGV0ZSBuYW1lIG9mIHRoZSBjb3VudHJpZXMuIApEaXNjIDwtIERpc2MgfD4gbXV0YXRlKGNvdW50cnkgPSBjYXNlX3doZW4oCiAgaXNvY250cnk9PSAiQkUiIH4gIkJlbGdpdW0iLAogIGlzb2NudHJ5PT0gIkRLIiB+ICJEZW5tYXJrIiwKICBpc29jbnRyeT09ICJHUiIgfiAiR3JlZWNlIiwKICBpc29jbnRyeT09ICJFUyIgfiAiU3BhaW4iLAogIGlzb2NudHJ5PT0gIkZJIiB+IkZpbmxhbmQiLAogIGlzb2NudHJ5PT0gIkZSIiB+ICJGcmFuY2UiLAogIGlzb2NudHJ5PT0gIklFIiB+ICJJcmVsYW5kIiwKICBpc29jbnRyeT09ICJJVCIgfiAiSXRhbHkiLAogIGlzb2NudHJ5PT0gIkxVIiB+ICJMdXhlbWJvdXJnIiwKICBpc29jbnRyeT09ICJOTCIgfiAiTmV0aGVybGFuZHMiLAogIGlzb2NudHJ5PT0gIkFUIiB+ICJBdXN0cmlhIiwKICBpc29jbnRyeT09ICJQVCIgfiAiUG9ydHVnYWwiLAogIGlzb2NudHJ5PT0gIlNFIiB+ICJTd2VkZW4iLAogIGlzb2NudHJ5PT0gIkRFLVciIH4gIkdlcm1hbnkgKHNwZWNpZmljYWxseSB0aGUgc3RhdGUgb2YgTm9ydGggUmhpbmUtV2VzdHBoYWxpYSkiLAogIGlzb2NudHJ5PT0gIkRFLUUiIH4gIkdlcm1hbnkgKHNwZWNpZmljYWxseSB0aGUgc3RhdGUgb2YgQmVybGluKSIsCiAgaXNvY250cnk9PSAiR0IiIH4gIlVuaXRlZCBLaW5nZG9tIiwKICBpc29jbnRyeT09ICJCRyIgfiAiQnVsZ2FyaWEiLAogIGlzb2NudHJ5PT0gIkNZIiB+ICJDeXBydXMiLAogIGlzb2NudHJ5PT0gIkNaIiB+ICJDemVjaCBSZXB1YmxpYyIsCiAgaXNvY250cnk9PSAiRUUiIH4gIkVzdG9uaWEiLAogIGlzb2NudHJ5PT0gIkhVIiB+ICJIdW5nYXJ5IiwKICBpc29jbnRyeT09ICJMViIgfiAiTGF0dmlhIiwKICBpc29jbnRyeT09ICJMVCIgfiAiTGl0aHVhbmlhIiwKICBpc29jbnRyeT09ICJNVCIgfiAiTWFsdGEiLAogIGlzb2NudHJ5PT0gIlBMIiB+ICJQb2xhbmQiLAogIGlzb2NudHJ5PT0gIlJPIiB+ICJSb21hbmlhIiwKICBpc29jbnRyeT09ICJTSyIgfiAiU2xvdmFraWEiLAogIGlzb2NudHJ5PT0gIlNJIiB+ICJTbG92ZW5pYSIsCiAgaXNvY250cnk9PSAiSFIiIH4gIkNyb2F0aWEiLAogIFRSVUUgfiBOQV9jaGFyYWN0ZXJfKSkKCnZpZXcoRGlzYykKCmBgYApDb252ZXJ0IG51bWVyaWMgYW5zd2VycyBmcm9tIHFjMTkgdG8gY2hhcmFjdGVyIGFuc3dlcnMKCmBgYHtyfQpEaXNjIDwtIERpc2MgJT4lIG11dGF0ZShjaGFuZ2VfZG9jcyA9IGNhc2Vfd2hlbigKICBxYzE5ID09IDEgfiAiWWVzIiwKICBxYzE5ID09IDIgfiAiTm8iLCAKICBxYzE5ID09IDMgfiAiREsiLCAKICBUUlVFIH4gTkFfY2hhcmFjdGVyXykpIApgYGAKCgpDb3VudCB0aGUgcWMxOSBhbnN3ZXJzIGJ5IGNvdW50cnkgCgpgYGB7cn0KbGlicmFyeShkcGx5cikKCiNzdW0gY291bnRzIApvcGluaW9uc19ieV9jb3VudHJ5IDwtIERpc2MgJT4lCiAgZ3JvdXBfYnkoY291bnRyeSwgY2hhbmdlX2RvY3MpICU+JQogIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lCiAgcGl2b3Rfd2lkZXIobmFtZXNfZnJvbSA9IGNoYW5nZV9kb2NzLCB2YWx1ZXNfZnJvbSA9IGNvdW50LCB2YWx1ZXNfZmlsbCA9IDApCgpvcGluaW9uc19ieV9jb3VudHJ5CgojcGVyY2VudGFnZSBjb3VudHMgCkRpc2MgJT4lCiAgZ3JvdXBfYnkoY291bnRyeSwgY2hhbmdlX2RvY3MpICU+JQogIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lCiAgdW5ncm91cCgpICU+JQogIGdyb3VwX2J5KGNvdW50cnkpICU+JQogIG11dGF0ZSh0b3RhbF9jb3VudCA9IHN1bShjb3VudCkpICU+JQogIG11dGF0ZShwZXJjZW50YWdlID0gKGNvdW50IC8gdG90YWxfY291bnQpICogMTAwKSAlPiUKICBzZWxlY3QoY291bnRyeSwgY2hhbmdlX2RvY3MsIHBlcmNlbnRhZ2UpICU+JQogIHBpdm90X3dpZGVyKG5hbWVzX2Zyb20gPSBjaGFuZ2VfZG9jcywgdmFsdWVzX2Zyb20gPSBwZXJjZW50YWdlLCB2YWx1ZXNfZmlsbCA9IDApCgpgYGAKYGBge3J9CiNUb3AgNSBvZiBjb3VudHJpZXMgdGhhdCB0aGluayB0aGF0IHRyYW5zZ2VuZGVyIHBlcnNvbnMgc2hvdWxkIGJlIGFibGUgdG8gY2hhbmdlIHRoZWlyIGNpdmlsIGRvY3VtZW50cyB0byBtYXRjaCB0aGVpciBpbm5lciBnZW5kZXIgaWRlbnRpdHkKY291bnRyeV93aXRoX21vc3RfeWVzIDwtIG9waW5pb25zX2J5X2NvdW50cnkgJT4lCiAgZmlsdGVyKFllcyA+IDApICU+JQogIGFycmFuZ2UoZGVzYyhZZXMpKSAlPiUKICBoZWFkKDUpCmNvdW50cnlfd2l0aF9tb3N0X3llcwpgYGAKCmBgYHtyfQojVG9wIDUgb2YgY291bnRyaWVzIHRoYXQgZG8gbm90IHRoaW5rIHRoYXQgdHJhbnNnZW5kZXIgcGVyc29ucyBzaG91bGQgYmUgYWJsZSB0byBjaGFuZ2UgdGhlaXIgY2l2aWwgZG9jdW1lbnRzIHRvIG1hdGNoIHRoZWlyIGlubmVyIGdlbmRlciBpZGVudGl0eQpjb3VudHJ5X3dpdGhfbW9zdF9ubyA8LSBvcGluaW9uc19ieV9jb3VudHJ5ICU+JQogIGZpbHRlcihObyA+IDApICU+JQogIGFycmFuZ2UoZGVzYyhObykpICU+JQogIGhlYWQoNSkKY291bnRyeV93aXRoX21vc3Rfbm8KYGBgCmBgYHtyfQojVG9wIDUgb2YgY291bnRyaWVzIHRoYXQgZG8gbm90IHRoaW5rIHRoYXQgdHJhbnNnZW5kZXIgcGVyc29ucyBzaG91bGQgYmUgYWJsZSB0byBjaGFuZ2UgdGhlaXIgY2l2aWwgZG9jdW1lbnRzIHRvIG1hdGNoIHRoZWlyIGlubmVyIGdlbmRlciBpZGVudGl0eQpjb3VudHJ5X3dpdGhfbW9zdF9ESyA8LSBvcGluaW9uc19ieV9jb3VudHJ5ICU+JQogIGZpbHRlcihESyA+IDApICU+JQogIGFycmFuZ2UoZGVzYyhESykpICU+JQogIGhlYWQoNSkKY291bnRyeV93aXRoX21vc3RfREsKYGBgCgo=